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We present a two stage hybrid inflationary scenario in non-minimal supergravity which can pre¬ 
dict values of the tensor-to-scalar ratio of the order of few x 10“^. For the parameters considered, 
the underlying supersymmetric particle physics model possesses two inflationary paths, the trivial 
and the semi-shifted one. The trivial path is stabilized by supergravity corrections and supports a 
first stage of inflation with a limited number of e-foldings. The tensor-to-scalar ratio can become 
appreciable while the value of the scalar spectral index remains acceptable as a result of the competi¬ 
tion between the relatively mild supergravity corrections and the strong radiative corrections to the 
inflationary potential. The additional number of e-foldings required for solving the puzzles of hot 
big bang cosmology are generated by a second stage of inflation taking place along the semi-shifted 
path. This is possible only because the semi-shifted path is almost perpendicular to the trivial one 
and, thus, not affected by the strong radiative corrections along the trivial path and also because 
the supergravity effects remain mild. The requirement that the running of the scalar spectral in¬ 
dex remains acceptable limits the possible values of the tensor-to-scalar ratio not to exceed about 
5 X 10“^. Our model predicts the formation of an unstable string-monopole network, which may 
lead to detectable gravity wave signatures in future space-based laser interferometer observations. 

PACS numbers: 98.80.Cq 


I. INTRODUCTION 

Inflation (for a review see e.g. Ref. Q) is by now 
considered to be an integral part of standard cosmology 
thanks to a plethora of precise observations on the cosmic 
microwave background radiation (CMBR) and the large- 
scale structure in the universe. Therefore, it is very im¬ 
portant to construct realistic inflationary models based 
on particle theory and consistent with all the available 
cosmological and phenomenological requirements. Un¬ 
doubtedly, hybrid inflation Q is one of the most promis¬ 
ing inflationary scenarios. It is ii naturally realized 
in the context of supersymmetric (SUSY) grand unified 
theory (GUT) models based on gauge groups with rank 
greater than or equal to five. 

In standard SUSY hybrid inflation, however, the GUT 
gauge symmetry is spontaneously broken only at the end 
of inflation and, thus, if magnetic monopoles are pre¬ 
dicted by this symmetry breaking, they are copiously 
produced @ , leading to a cosmolc^ical catastrOThe. This 
disaster is avoided in the smooth or shifted [g variants 
of SUSY hybrid inflation, where the GUT gauge symme¬ 
try is broken already during inflation. These variants 
were based on non-renormalizable superpotential terms. 
It was, though, subsequently shown that a new smooth 
0 and a new shifted Q hybrid inflation scenario can 
be constructed with only renormalizable superpotential 
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terms within an extended Pati-Salam (PS) SUSY GUT 
model, which was initially introduced Q for solving a 
very different problem. Namely, the simplest SUSY PS 
model predicts (see Ref. 0) exact Yukawa unification 
m and, if it is supplemented with universal bound¬ 
ary conditions, yields unacceptable 6-quark mass values. 
This problem is solved in the extended model, where 
Yukawa unification is naturally and moderately violated. 

After the first accurate measurement 0 of the scalar 
spectral index rig, however, it has been realized that there 
is a tension between all these well-motivated and natu¬ 
ral inflationary scenarios and the measured value of this 
index. Indeed, within the standard power-law cosmo¬ 
logical model with cold dark matter and a cosmological 
constant, the data imply that Us is clearly lower than 
unity - for the latest results on rig see Ref. [1^ . Inflation¬ 
ary scenarios, on the other hand, such as the ones men¬ 
tioned above, within super grav ity (SUGRA) with mini¬ 
mal Kahler potential, yield [1^ rig’s which are very close 
to unity or even exceed it. 

One idea [l^ for reducing the predicted spectral index 
is based on the observation that rig generally decreases 
with the number of e-foldings suffered by our present 
horizon scale during inflation. So, reducing this num¬ 
ber of e-foldings, we can achieve values of rig compatible 
with the recent data without having to abandon the use 
of a minimal Kahler potential. The additional number 
of e-foldings required for solving the horizon and flat¬ 
ness problems of standard hot big bang cosmology can 
be provided by a subsequent second stage of inflation. 

It is interesting to note that the extended SUSY PS 
model of Ref. Q, which can lead to new smooth Q or 
new shifted @ hybrid inflation, can also provide us with a 
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double inflation scenario called standard-smooth hybrid 
inflation [T^ which solves the above mentioned spectral 
index problem along the lines just discussed. The cos¬ 
mological scales exit the horizon during the main stage 
of inflation, which is of the standard hybrid type and 
occurs as the system slowly rolls down a trivial classi¬ 
cally flat direction on which the PS gauge group is un¬ 
broken. This direction is subsequently destabilized giv¬ 
ing its place to a classically non-flat valley of minima 
along which new smooth hybrid inflation takes place with 
the PS GUT gauge group being broken. Consequently, 
magnetic monopoles are produced only at the end of the 
first stage of inflation, but they are adequately diluted by 
the second stage, which also provides the extra e-foldings 
needed for solving the puzzles of hot big bang cosmology. 

After the recent results of BICEP2 [a on the B-mode 
in the polarization of the CMBR at degree angular scales, 
it seems possible that the inflationary scenarios will have 
to face a new challenge. Namely, they should be able to 
accommodate appreciable values of the tensor-to-scalar 
ratio r, since a B-mode of primordial origin could be 
due to the production of gravitational waves during infla¬ 
tion. We should, however, consider this possibility with 
reservation since some serious criticism [l8| to the origi¬ 
nal BICEP2 analysis has already appeared claiming that 
the foreground from Galactic polarized-dust emission has 
been underestimated. On the other hand, after the re- 
centl y re leased Planck HFI 353 GHz dust polarization 
data [^, the first attempts to make a joint analysis of the 
Planck and BICEP2 data have been presented [10, |2l|. 
They showed that, although r is smaller than initially 
claimed, significant values of r - of order 0.01 - cannot 
be excluded. The most recent joint analysis [Uj yields an 
upper limit on r of about 0.12 at 95% confidence level. 
Unfortunately, all the above mentioned variants of SUSY 
hybrid inflation predict negligible values of r. So, it is cer¬ 
tainly worth investigating whether realistic SUSY hybrid 
inflation models accommodating appreciable values of r 
can be constructed. 

In Ref. [13, a double inflation scenario has been pro¬ 
posed which is compatible with the BICEP2 data [l7l |. 
The first stage of inflation is of the SUSY hybrid type, 
while the second stage is left unspecified, which makes 
the scenario incomplete. The inflationary potential is 
supplemented with a mass-squared term for the inflaton 
attributed to SUGRA corrections and with a logarith¬ 
mic term representing very strong radiative corrections 
due to the SUSY breaking during inflation. It is the 
competition between these two contributions which al¬ 
lows appreciable values of r while Ug remains acceptable. 
The assumption, however, that an inflaton mass-squared 
term is the only relevant SUGRA correction during infla¬ 
tion in a scenario with Planck-scale inflaton field values 
seems totally unjustified. In addition, this paper follows 
the usual practice of only taking into account the ra¬ 
diative corrections in the derivatives of the inflationary 
potential when calculating the slow-roll parameters and 
neglecting them when calculating the potential itself. In 


the case of extremely strong radiative corrections, how¬ 
ever, this may lead to erroneous results. Also, Ref. [^ 
attempts to accommodate the BIGEP2 results in a dou¬ 
ble hybrid inflation model where the inflaton potential 
changes dynamically with the evolution of the inflaton 
fields. The particular implementation of this interesting 
idea, though, appears to have some problems of natural¬ 
ness in the design of the superpotential. In addition, the 
treatment of SUGRA seems to be incomplete. Finally, 
Ref. 2M shows that, in SUSY hybrid inflation models, it 
is possible to obtain values of r close to 0.03 by employing 
an expansion of a non-minimal Kahler potential with ap¬ 
propriate coefficients. The validity of this approach may, 
however, be questionable since the inflaton takes values 
close to the Planck scale. 

In this paper, we will show that a reduced version of the 
extended SUSY PS model of Ref. (3 can yield a two stage 
inflationary scenario which can predict values of r up to 
about 0.05. Larger values of the tensor-to-scalar ratio 
would lead to unacceptably large running of the scalar 
spectral index. In the range of the model parameters 
considered here, the model in the global SUSY limit pos¬ 
sesses practically two classically flat directions, namely 
the trivial and the semi-shifted [13 one. After including 
SUGRA corrections, the trivial path, on which the full 
GUT gauge group is unbroken, is stabilized and a first 
stage of inflation can occur as the system slowly rolls 
down this path. All the cosmological scales exit the hori¬ 
zon during this stage and our present horizon undergoes 
a limited number of e-foldings. The obtained tensor-to- 
scalar ratio can be appreciable while the scalar spectral 
index assumes acceptable values thanks to the compet¬ 
ing effect of the sufficiently mild SUGRA corrections re¬ 
sulting from the construction of Ref. [l^] and the strong 
radiative corrections to the inflationary potential. 

Subsequently, a second inflationary stage occurs along 
the semi-shifted path, where U(1)b-l remains unbroken, 
and provides the additional number of e-foldings required 
for solving the standard problems of hot big bang cosmol¬ 
ogy. This is possible since, for our choice of parameters, 
the semi-shifted path is almost perpendicular to the triv¬ 
ial path and, thus, is not affected by the strong radiative 
corrections on the trivial path. It is also important that 
the SUGRA corrections on the semi-shifted path are kept 
sufficiently mild again by the mechanism of Ref. [l^l. 

We first present, in Sec. m the salient features of the 
model in global SUSY. In Sec. IIIIl we then calculate the 
SUGRA and one-loop radiative corrections to the poten¬ 
tial and discuss our double inflationary scenario. Finally, 
in Sec. CYl we summarize our conclusions. Through¬ 
out, we will use units where the reduced Planck mass 
Top = 2.4354 x 10^® GeV is equal to unity. 


II. THE MODEL IN GLOBAL SUSY 

We consider a reduced version of the extended SUSY 
PS model of Ref. Q- This version is based on the left- 
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right symmetric gauge group Glr = SU{3)c x S'C/(2)l x 
SU{2)^xU{1)b-l, which is a subgroup of the PS group. 
The superfields of the model which are relevant for in¬ 
flation are the following. A gauge singlet S, a pair of 
superfields $, l> belonging to the ( 1 , 1 ,3)o representa¬ 
tion of Glr, and a conjugate pair of Higgs superfields 
H and H belonging to the (1,1, 2)i and (l,l,2)_i rep¬ 
resentations of Glr, respectively. The field $ acquires 
a vacuum expectation value (VEV) which breaks Glr 
to Gsm X U{1)b-l with Gsm being the standard model 
(SM) gauge group, while the VEVs of H and H cause 
the breaking of Glr to Gsm- The full superfield con¬ 
tent and superpotential, the global symmetries, and the 
charge assignments can be easily derived from the ex¬ 
tended SUSY PS model of Ref. Q by simply reducing its 
GUT gauge group to Glr. The only global symmetry of 
the model which is relevant here is its U{1) R symmetry 
under which S and $ have charge 1 with all the other 
superfields mentioned above being neutral. 

The superpotential terms relevant for inflation are 

W = kS {M'^ - --/SHH+ m^^-X^HH, (1) 

where M, m are superheavy masses and k, 7 , A are di¬ 
mensionless coupling constants. These parameters are 
normalized so that they correspond to the couplings be¬ 
tween the SM singlet components of the superfields. The 
mass parameters M, m and any two of the three di¬ 
mensionless parameters k, 7 , A can always be made real 
and positive by appropriately redefining the phases of 
the superfields. The third dimensionless parameter, how¬ 
ever, remains generally complex. For definiteness, we will 
choose this parameter to be real and positive too. 

The F-term scalar potential obtained from the super¬ 
potential in Eq. o is given by 

= |k(M2 - $2) _ 

+ \m^ - 2kS'$|^ -b |to$ - XHH\'^ 

+ |7,5 + Ai|2(|i/p + |i7|2), (2) 

where the complex scalar fields which belong to the SM 
singlet components of the superfields are denoted by the 
same symbol. From this potential and the vanishing of 
the D-terms (which implies that H* = one finds 

[ 2 ^ two distinct continua of SUSY vacua: 


with 


V° = Utr = (7) 


On this direction, Glr is unbroken. The second one, 
which appears at 


$ = - 
HH = 

= Ksh = 


jm 


^ s 

A ’ 


2kX’ 

— d)^) + Amd? 


7" 


A 2 


A 2 


1 -b 


2 2 \ 
7 771 \ 


4k2A2M2 


1 , ( 8 ) 


is the trajectory for the new shifted hybrid inflation Q. 
On this direction, Glr is broken to Gsm- The third flat 
direction, which exists only if > rn?/2 k^, lies at 


d- = ±M 



m? 

2 k 2 M 2 ’ 


d> = 


— 5. 

m 


H = H = 0. 


(9) 

It is the path along which semi-shifted hybrid inflation 
[ 2 ^ takes place with 


Ui> = U,,h = 


( rn? 
V ~ 4k2M2 


( 10 ) 


Along this direction Glr is broken to Gsm x U{\)b-l- 
We choose to consider the case where > rv?I2i^ 
and, thus, the semi-shifted flat direction exists. One can 
show - see Ref. [ 2 ^ - that, in this case, we always have 
Ussh < Uish and Ussh < Utr- Therefore, the semi-shifted 
flat direction, if it exists, always lies lower than both the 
trivial and the new shifted one. On the other hand, the 
new shifted flat direction may either lie lower or higher 
than the trivial one depending on the values of the pa¬ 
rameters. Here we will take 7 ,cA<Ck, m<C M, 

and [S'! < 1. In this case, the new shifted flat direction 
practically coincides with the trivial one and, thus, plays 
no independent role in our scheme. 


III. THE DOUBLE INFLATIONARY SCENARIO 


$ = d>+, H* = H, \H\ = 5 = !> = 0, (3) 

$ = H* = -H, \H\ = 5 = $ = 0, (4) 

where 

The potential in Eq. generally, possesses [ 1 ^ three 
flat directions. The first one is the usual trivial flat di¬ 
rection at 

$ = |. = 77 = iJ = 0 ( 6 ) 


In this section, we will show that, after including 
SUGRA corrections, the trivial path becomes stable for 
large absolute values of the real canonically normalized 
inflaton. Thus, it can support a first stage of infla¬ 
tion during which the universe undergoes a number of 
e-foldings which, although limited, is adequately large 
for all the cosmological scales to exit the horizon. Strong 
radiative corrections to the inflationary potential, which 
are controlled by the parameter k, in conjunction with 
mild SURGA corrections then guarantee that an appre¬ 
ciable value of the tensor-to-scalar ratio can be achieved 
together with an acceptable value of the scalar spectral 
index. Our scenario can predict values of the tensor-to- 
scalar ratio only up to about 0.05 because larger values 
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require unacceptably large running of the scalar spectral 
index. 

A subsequent second stage of inflation along the semi- 
shifted path can provide the additional number of e- 
foldings required for solving the horizon and flatness 
problems of the standard hot big bang cosmology. This 
is possible since, for the parameters chosen, this direction 
is almost orthogonal to the trivial path and, thus, it is 
not affected by the strong radiative corrections present 
during the first stage of inflation. In this connection, 
it is also important that the SUGRA corrections on the 
semi-shifted path remain mild. 

After the termination of the first stage of inflation, 
the system moves towards the semi-shifted path and 
the group SU{2)b. breaks spontaneously to a ?7(1) sub¬ 
group, leading to the formation of magnetic monopoles. 
On the other hand, the spontaneous breaking of a lin¬ 
ear combination of this C/(l) and U{1)b-l, which takes 
place at the end of the second inflationary stage, leads 
to the production of open cosmic strings which connect 
these monopoles to antimonopoles. Subsequently, the 
monopoles come into the post-inflationary horizon and 
the whole system of strings and monopoles decays well 
before recombination without leaving any trace in the 
CMBR. The gravitational waves which are generated by 
the decaying strings may, though, be measurable in the 
future. 


A. The first inflationary stage 


parametrized by the complex inflatons S and $_( approx¬ 
imately). However, as we shall see - cf. Ref. - , the 
relation 

\Kz,\^K^^z,^ = \Z2\^ = P (14) 

implies that these inflatons acquire squared masses pro¬ 
portional to P as soon as the value of Z 2 becomes non¬ 
zero. 

Using the R and U{1)b-l symmetries of the model, we 
can rotate S and H on the real axis - cf. e.g. Ref. [^ . 
The fields $, <I>, H remain in general complex. However, 
for simplicity, we will also restrict them on the real axis. 
This is not expected to influence our results in any essen¬ 
tial way since these fields are anyway real in the vacuum 
and on all the flat directions of the model given that 
the parameters of the model are chosen to be real - see 
Sec. EH Also, we can show that, everywhere on the trivial 
and the semi-shifted inflationary paths, the mass squared 
matrices of the imaginary parts of the fields do not mix 
with the mass squared matrices of their real parts and, 
during both inflations, have positive eigenvalues in the 
directions perpendicular to these paths. So there is no 
instability in the direction of the imaginary parts of the 
fields which are orthogonal to these inflationary trajecto¬ 
ries. Moreover, as we can demonstrate, both the trivial 
and the semi-shifted inflationary paths are destabilized 
with the fields developing real values. 

We define the canonically normalized real scalar fields 
cr, (j), (j), h, h corresponding to the Kahler potential in 
Eq. (fTTIl as follows - cf. Ref. [ 2 ^ - : 


We adopt here the following Kahler potential 

K = - In (1 - |5|2) - In (1 - |l>|2) + + \H\^ 

+ |iJ|2-21n(-ln|Zip) + |Z2p (11) 

{\s\, |l>| < 1, 0 < \Zi\ < 1), where we included two extra 
Glr singlet superfields Zi and Z 2 , which do not enter the 
superpotential at all because they transform non-trivially 
under additional anomalous 1/(1) gauge symmetries. The 
resulting F-term potential in SUGRA is given by 


Vf 




3|1U|^ 


e^, (12) 


where a subscript denotes derivation with respect to 
the field Xi and the sum extends over all the seven fields 
S', l>, $, iJ, H, Zi^ 2 - The values of Z\ and Z 2 are as¬ 
sumed to be fixed by anomalous D-terms. Note that 
the superhelds S, l>, Z\ possess Kahler potentials of the 
no-scale type which for Z 2 = 0 , in view of the relation 

=2, (13) 

guarantee the exact flatness of the potential along the 
trivial path and its approximate flatness on the semi- 
shifted one - see below. These paths are, respectively. 


S = tanh —$ = tanh — 

72 72 


(15) 





(16) 


We can now evaluate the potential Vp in Eq. m 
with the overall factor exp [—2 In (—In |Zip) -|- IZ 2 P] ab¬ 
sorbed into redefined parameters k, 7 , m, and A and find 


Vf = 


A\ cosh^ — A 2 sinh^ -|- /SAg -|- A 4 -|- A 5 

+ + V) Al + \(4,^ + h-‘ + V)Al 

-h (72(/)A5 - 2hhA(^ Ag (17) 


Here 


Ai = K ( — 2— 1 — 'xhh. 


(18) 


A ^ ^hh 

A 2 = m—= - fifi 

72 2 


(19) 


Ag = Ai sinh cosh -|- A 2 cosh sinh , (20) 
v2 v2 v2 y2 
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A 4 = Ai sinh sinh —p + A 2 cosh cosh —p, ( 21 ) 

v2 \/2 y2 y2 


Ac. = m cosh —= sinh —■= — \/ 2 K(b sinh —■= cosh —=, 

\/2 V 2 V 2 


and 


( 22 ) 


Ag = 7 sinh cosh —p + A cosh sinh —=. (23) 

v2 v2 y2 v2 

On the trivial trajectory where (j), h, h = 0, the F- 
term potential takes the form 


Vf = 


1 + /3 sinh^ 


72. 


(24) 


The mass-squared eigenvalues in the directions perpen¬ 
dicular to this trajectory for sinh^ (ct/-\/ 2 ) M ^/2 can 

also be found from Eq. (113 to be 


~ sinh^ —p, 

^ ^/2 


(25) 


~ ( 1 -I- (1 -I- /3) sinh^ ) , (26) 


72 


m 


2 

XI 


and 




+ ((1 -I- I3)kM^ — 7 ) sinh^ 


72 J’ 

(27) 


-I- ((1 -I- I3)kM^ + 7 ) sinh^ , 

72 . 
(28) 

where Xi = {h + h)ly/2 and X 2 = {h — h)ly/2. Thus, 
assuming that 7 < we see that the trivial path, 

which is flat in the limit /3 —>■ 0, is stable for large absolute 
values of the inflaton a. 

Note that Eqs. dUl) and (1^51) hold for any value of 
sinh^ ( 17 / 72 ). On the contrary, one can show that as 
sinh^ (tT/ 72 ) decreases, the eigenvalues and eigenstates 
of the mass-squared matrix of the (j)—4> system change. In 
particular, when sinh^ ((t/72) ~ M^/2 -|- m?/{ 2 k^M'^), 
the mass-squared matrix of the </ — </> system acquires 
a zero eigenvalue with </> dominating the corresponding 
eigenstate. Subsequently, as sinh^((T/v^ approaches 
the value M^/2, the eigenvalues of the (/—</ mass-squared 
matrix become almost opposite to each other with (j) and 
(j) contributing almost equally to both the eigenstates. 
A further decrease of sinh^ (cr/7^ leads to the domina¬ 
tion of the unstable eigenstate by </>. Since the field (j) 
is required to develop a nonzero VEV in order to can¬ 
cel the false vacuum energy density k^M'^ on the trivial 
trajectory - see Eq. @ or Eqs. (HB and m- , we will 




take as critical value CTc of a at which the trivial path is 
destabilized the one determined by the relation 


sinh^ 

72 


M 2 


(29) 


To the F-term scalar potential Vp in Eq. (II3 has to 
be added during the first stage of inflation (i.e. for \a\ > 
|(7c|) the term 


vf = 


f 


In ■ 


2 tanh^ 


V 2 


M 2 


(30) 


corresponding to the dominant one-loop radiative cor¬ 
rections to the inflationary potential due to the N^- 
dimensional supermultiplet $ (fV^ = 3). Notice that the 
renormalization scale in these radiative corrections is cho¬ 
sen such that vf vanishes at |cr| = IcTcl (tanh^ (ctc/T^ — 
sinh^ (ctc/T^ = M 2 / 2 ). 

Setting 


54 , = 


N^k'^ 

27r^ 


(31) 


we can rewrite the full inflationary potential and its 
derivatives (denoted by primes) with respect to the 
canonically normalized real inflaton field a as follows: 


. =14-/3 sinh^ ^ In 
k2M4 a/2 4 


2 tanh' 


2 (7 


V' 

k‘^M^ 


= sinh(72CT) (/3 -I- 

72 V si 


M 2 


71 ^ da), 


sinh^(72cr) 


(32) 

(33) 


V'‘ 


= cosh(72cr) (/3-^^ 


and 


V" 


2M4 


sinh^(72<7) 


= 72sinh(72(7) (/3-;— 


(34) 


sin7(72cr) 


-4 


‘2^/254 


tanh^(72CT) sinh(72i7) 


(35) 


The usual slow-roll parameters for inflation are then writ¬ 
ten as 


e = 


V' ^ 


(36) 

j 

I C^(a)’ 

V" ^ 

) 1 

(37) 

k2M4^ 

> c{<yy 


and 
? = 


«2M4 ) ) (72 (o-) 

_4(30 y/e_ 

C{a) tanh^(72cr) |sinh(72cr)| 


tanh(72cr) t/Tc 


( 38 ) 
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Using these expressions, we can evaluate the scalar spec¬ 
tral index rig, its running Og, and the tensor-to-scalar 
ratio r from the formulas 

rig = 1 -I- 2ry — 6 e, (39) 

Og = 167;e - 24e2 - 2^, (40) 

r = 16e. (41) 

Finally, the scalar potential on the trivial inflationary 
path can be written in terms of the scalar power spectrum 
amplitude Ag and r as follows 

U=—Agr. (42) 

As a numerical example, we take the value of the real 
inflaton field cr at horizon exit of the pivot scale fc* = 
0.05 Mpc“^ to be cr* = 1.45. Also, we take k = 1.7, 
P = 0 . 022 , and the scalar power spectrum amplitude 
Ag = 2.215 X 10“® at the same pivot scale [l^. With 
these input numbers, we then find M = 3.493 x 10“^, 
C{a^) = 2.2941, e = 0.00188, rj = -0.01389, rig = 0.9609, 
Og = —0.01674, and r = 0.0301. 

As one can see, our predictions can not only be per¬ 
fectly consistent with the latest data released by the 
Planck satellite experiment [I^, but also accommodate 
large values of the tensor-to-scalar ratio r of order few x 
10“^. As it is obvious from Eq. dUD, such values of r re¬ 
quire relatively large values of e, which in turn reduce the 
scalar spectral index rig in Eq. (1391) below unity, but not 
quite adequately to make it compatible with the data. So 
we need an appreciable negative value of rj, which requires 
that the parenthesis in the right hand side of Eq. (IMl) be 
dominated by the second term. A similar parenthesis ap¬ 
pears in the right hand side of Eq. (1551) too, but with the 
two terms added rather than subtracted. As it turns out, 
both these terms have to be appreciable with the second 
one being larger in order to be able to bring rig near its 
best-fit value from the Planck data. This is certainly pos¬ 
sible only for large values of the parameter k controlling 
the radiative corrections on the trivial path. Note that 
the first stage of inflation ends before the system reaches 
the critical point in Eq. (1551) by violating the slow-roll 
conditions and the obtained number of e-foldings is lim¬ 
ited due to the large values of e involved and the fact 
that CT* ^ 1 . 

B. The second inflationary stage 

For the rest of the parameters of the model, we chose 
the values m = 1.827 x 10“^, A = 0.1, and 7 = 10“®. We 
solved numerically the differential equations of the sys¬ 
tem with potential energy density given by the exact Vp 
in Eq. m supplemented with the relevant radiative cor¬ 
rections and the D-terms involving the fields H, H. The 


numerical investigation then revealed that there exist ap¬ 
propriate small initial absolute values of the scalar fields 
if), (j), h, h for which, after the first stage of inflation and 
the elapse of a sufficient amount of cosmic time t for the 
energy density to approach rn?M‘^, we have (jP ~ 2 M^, 
h, h ~ 0 , and the scalar fields cr and 0 take values such 
that A 5 ~ 0 with |cr| <C 1. So it is obvious that the sys¬ 
tem reaches the semi-shifted inflationary path in Eq. m 
- note that the second relation in this equation is equiva¬ 
lent to A 5 = 0. It is remarkable that |0|, which at the end 
of the first inflationary stage is extremely small, manages 
to attain values of the order of few x 10 “^ at the onset 
of the second stage. 

It is worth noticing, however, that the initial values of 
the fields which lead to a double inflation scenario, al¬ 
though quite frequent, do not seem to form well-defined 
connected regions. In other words, the solutions of the 
coupled system of differential equations exhibit a rather 
‘chaotic’ behavior in the sense that a slight change of 
the initial conditions may lead from a double to a sin¬ 
gle inflation scenario. Note that a similar situation is 
encountered even in the simplest SUSY hybrid infla¬ 
tion scenario, where a slight change of initial conditions 
may ruin inflation leading the system directly to the vac¬ 
uum. Such a behavior, given the multidimensionality of 
the field space in our case, makes it very difficult to pro¬ 
vide further details concerning the structure of the space 
of initial values leading to the desirable scenario. 

For a negligible value of 7 , (jP ~ 2M^, A 5 ~ 0, and 
|cr| <C I, we find that Ai ~ 0 , A 3 ~ A 2 sinh ( 0 /-\/ 2 ), 
A 4 ~ A 2 cosh ((^/v^), Ae ~ Asinh (0/-\/2), and the F- 
term scalar potential becomes 


Uf 


^2 + (/3 + M^)A2 sinh^ 

v2 

+ -{hP + JP) (A^ -I- A2) smh^ 


— 2 / 1 / 1 AA 2 sinh^ 


V2J 


M^/3 




1 + P sinh^ 


(44) 


72. 


The expression in Eq. (1441) gives approximately the E- 
term potential on the semi-shifted path. Notice the strik¬ 
ing similarity of this expression with the one in Eq. (1241) 
involving the same parameter p. 

From A 5 ~ 0 and the fact that A 5 oc mtanh (0/72) — 
72K0tanh (cr/7^, it follows that the combination of S 
and $ which could remain large when the energy den¬ 
sity approaches irpM"^ and plays the role of the complex 
inflaton in the second stage of inflation is 


mS -I- 2k < <1) > 4) ^ 

'Pm? 4k^M^ 


since the contribution of $ in this combination is about 
2KM/m ~ 650 times bigger than the one of S. 
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From Eq. (I43|) , we can construct the mass-squared ma¬ 
trix for the h — h system during the second stage of in¬ 
flation. We find that the mass eigenstates are given by 
the combinations Xi = (/i -I- h)/'/2 and X 2 = {h — h)l\/2 
with masses-squared 

(A — (1 -I- (3)mM) sinh^ — mM , 
v 2 J 
(46) 

(A -I- (1 -I- (3)mM) sinh^ + mM . 

v 2 

.(47) 

We see that xi develops an instability which terminates 
the valley along which the second stage of inflation takes 
place with the critical value 4>c of the real canonically 
normalized inflaton (p being approximately determined 
from the relation 


— (-^ + mM) 


~ (A — mM) 


sinh^ 


0 c 

71 


mM 


(48) 


have seen, for accommodating appreciable values of r, do 
not affect the second stage of inflation. This is very cru¬ 
cial since otherwise the semi-shifted path would become 
too steep and there would be no way of generating the 
extra e-foldings required for solving the puzzles of hot big 
bang cosmology. 

The number of e-foldings during the second stage of 
inflation between an initial value 0 in and a final value 0 f 
of the inflaton 0 is given, in the slow-roll approximation, 
by iV(0f) - iV(0in), where 


N(4,) = ‘ In + y /l - (j^ 

2/3V1 - (4/0) cosh(^0) - ^1 - (4/0) 

(51) 

with 


^ NhX^ 


(52) 


The termination of slow-roll inflation is due to the radia¬ 
tive corrections in Eq. o and takes place at a value 0 f 
(|0f| > |0c|) of 0 given by 


To the F-term scalar potential Vp during the second 
stage of inflation (i.e. for | 0 | > | 0 c| and |cr| < |ctc| ) has 
to be added the following term 


yh. ^ 


m M ' 




(tanh -f tanh 

X In ('49') 


(l + 4«^S)m 


which may be approximated as 




^M- 


(NhX'^X 
\ J 


In- 


A tanh^ 


^/2 


mM 


(50) 


and corresponds to the dominant one-loop radiative cor¬ 
rections due to the TV/j-dimensional supermultiplets H, H 
{Nh = 2). Notice that the renormalization scale is cho¬ 
sen such that vanishes at |0| = |0c| (tanh^ (0c/'\/2) ~ 
sinh^ (0c/-\/2) = mM/X). 

The one-loop radiative corrections involving the $ su- 
permultiplet are neglected since they are relatively very 
small. This is because $ couples to the combination 
which plays the role of the complex inflaton during the 
second stage of inflation only through S and the contri¬ 
bution of S to this combination is severely suppressed. 
Indeed, the slope of the potential along the semi-shifted 
path generated by the radiative corrections involving the 
<i> supermultiplet is suppressed relative to the one involv¬ 
ing the H, H supermultiplets by, approximately, a factor 
(N^/SNh) {m/XM)^ ^ 5x10“"^. This is a very important 
property of our model resulting from the fact that, for the 
parameters chosen, the semi-shifted path is almost per¬ 
pendicular to the trivial one. So the very strong radiative 
corrections on the trivial trajectory, which are controlled 
by the strong coupling constant k and are needed, as we 


cosh(-/ 20 f) - Y + + Y' (^^) 

It turns out numerically that, with the chosen values of 
the parameters, the pivot scale /c* = 0.05 Mpc“^ suffers 
about 13 e-foldings during the first stage of inflation. As 
a consequence, approximately another 38-39 e-foldings 
(for reheat temperature = 10® GeV) must be pro¬ 
vided by the second inflationary stage, which requires a 
value of |0in| — 0.23 at the onset of this stage. This 
requirement can indeed be fulfilled in our numerical ex¬ 
ample as we have shown by extensive numerical studies. 
It is worth noticing that, due to the presence of mild but 
appreciable SUGRA corrections and not too weak radia¬ 
tive corrections, the second stage of inflation is able to 
generate a relatively limited number of e-foldings. Gon- 
sequently, this number is not too sensitive to the value 
of 0i„. 

In order to support the statements made above con¬ 
cerning the numerical part of our work, we depict, in 
Fig. E the evolution of the fields a and 0 as functions 
of the number of e-foldings N starting from the point 
where the pivot scale fc, exits the horizon for a particu¬ 
lar choice of initial conditions at this point. Namely, we 
start with a = 1.45, 0 = 10“^, 0 = 10“®, h = 10“^, and 
/i = 1.01 X 10“'’' corresponding to an almost D-flat direc¬ 
tion. All the fields are given zero initial velocity except 
for a the velocity of which is taken to be —1.1074 x 10“®. 
This is, indeed, the actual value of the velocity of cr on 
the trivial path, which was determined numerically. 

We observe that a assumes values above its critical 
value for about 13 e-foldings. Towards the end of the 
first inflationary stage, cr oscillates with appreciable am¬ 
plitude around zero four times. We may consider that, 
during these oscillations, the first stage of inflation has 
not come to an end. When the amplitude of the os¬ 
cillations falls below the critical value of cr, 0 moves to 



















FIG. 1 : The evolution of the fields a and (j> for the case with 
r — 0.0301 versus the number N of e-foldings after the horizon 
exit of the pivot scale fc», where we take a — 1 . 45 , 4 > = 10 ~®, 
0 = lO"®, h = 10 “'‘, h = 1.01 X 10 "'‘, and dajdt = - 1.1074 x 
10 "®. 



FIG. 2 : The evolution of the fields a and (j> for the case with 
r = 0.0502 versus the number N of e-foldings after the horizon 
exit of the pivot scale A:*, where we take a = 1 . 35 , (j) = 10 ~®, 
(f) = 10 ~®, h = 9 X 10 ““^, h = 9.01 X 10 ~^, and da/dt = 
- 1.8523 X 10 "®. 


its value on the semi-shifted path and (j) starts perform¬ 
ing slow oscillations with variable amplitudes typically of 
order M. The size of (j) remains small for about 1.7 e- 
foldings before starting its growth and acquires its largest 
value ~ 0.225 at iV ~ 17.7, when the second inflationary 
stage may be considered as having already started. For 
N > 20, the evolution of ^ follows Eq. (ED closely. 

If we allow for a stronger running of the scalar spectral 
index, we may obtain larger values of r. For example, 
taking cr* = 1.35, k = 1.75, and /3 = 0.037, we find 
M = 3.891 X 10-3, C(a^) = 2.3479, e = 0.00314, rj = 
-0.00844, Us = 0.9643, = -0.03007, and r = 0.0502. 

In addition, we choose m = 3.891 x lO”®, A = 0.1, and 
y = 10“®. With these choices, the pivot scale suffers 
about 10 e-foldings during the first stage of inflation and, 
consequently, approximately another 41 — 42 e-foldings 
must be provided by the second stage. This implies that 
l^inl lies in the range 0.38 — 0.40. We verified numerically 
that the fulfillment of this requirement is indeed feasible. 
In Fig. [2 we depict the evolution of the fields a and ^ 
as functions of the number of e-foldings N again starting 
from the point where the pivot scale fc, exits the horizon 
for a particular choice of initial conditions at this point. 

Needless to say that, changing the values of the input 
parameters of the model, we can easily achieve success¬ 
ful solutions with smaller values of the tensor-to-scalar 
ratio and that there is no particular fine-tuning of the 
parameters required in our model. 


C. The formation of monopoles and cosmic strings 

After the termination of the first stage of inflation, 
the system moves towards the semi-shifted path and 


the group S'C/(2 )r breaks spontaneously to a 17(1) sub¬ 
group by the nonzero value which the field 4> devel¬ 
ops. This leads to the formation of superheavy magnetic 
monopoles. We can obtain a rough order of magnitude 
estimate of the mean distance between the monopoles 
and the antimonopoles, which is though adequate for 
our purposes here. We assume that, at production, 
this distance is p(2kM)~^ as determined by the rele¬ 
vant Higgs boson mass with p ^ 1 being a geometric 
factor. In the matter dominated period between the 
two inflationary stages, this distance increases by a fac¬ 
tor ^ where k^M'^ and rri^M^ are 

the classical potential energy densities on the trivial and 
the semi-shifted paths, respectively. The subsequent sec¬ 
ond inflationary stage stretches this distance by a factor 
exp A^ 2 , where N 2 is the number of e-foldings du ring this 
stage which is large but not huge - cf. Ref. 28|. In 
the period of damped oscillations of the inflaton field, 
the monopole-antimonopole distance increases by an¬ 
other factor ^ /c{T,;)T^y/^, where Tj- is the re¬ 

heat temperature, which we take to be 10® GeV, and 
c(T) = 7r®5(T)/30 with g{T) being the effective number 
of massless degrees of freedom at cosmic temperature T. 
In the radiation dominated period which follows reheat¬ 
ing, the monopole-antimonopole distance is multiplied 
by another factor ^ T-^/T ^ (4c(T)/3)^/‘^Ti.-\/t, where t 
is the cosmic time. So this distance, in the radiation 
dominated period, becomes 



c{T,)-ic{T)ip 


(^) 


1 


T^t^. 


(54) 

Equating this distance with the post-inflationary parti¬ 
cle horizon ^ 2t, we find the time fn (in ip = ~ 
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2.7027 X 10 sec units) at which the monopoles enter 
this horizon: 


in 



c{Tr) 


^c(rH) 





(55) 


where Tjj is the cosmic temperature at cosmic time tn- 
The formation of monopoles is not the whole story 
though since our scenario leads to the generation of cos¬ 
mic strings too. Indeed, after the end of the second in¬ 
flationary stage, the system settles in one of the two dis¬ 
tinct continua of SUSY vacua in Eqs. and O with 
~ ±M in our case and a linear combination of the 
U(1)b-l gauge symmetry and the unbroken U(l) sub¬ 
group of SU{2 )b breaks spontaneously leading to the 
production of local cosmic strings. These strings, if they 
survived after recombination, could have a small contri¬ 
bution to the CMBR power spectrum parametrized [2^ 
by the dimensionless string tension Gfig, where G is New¬ 
ton’s gravitational constant and jig is the string tension, 
i.e. the energy per unit length of the string. Apply¬ 
ing to our case the results of Ref. [2^ , which considered 
local strings within the Abelian Higgs model in the Bo- 
gomol’nyi limit, we write, for the string tension, 

^ig = 4^1 (H) |2, (56) 


where (H) is the VEV of H. Although the strings in 
our model are more complicated than in the Bogomol’nyi 
limit of the Abelian Higgs model, we think that the above 
estimate for the string tension is good enough for our 
purposes here. 

In our case, the strings decay well before recombination 
and, thus, do not affect the CMBR. The reason is that 
they are mostly open strings connecting monopoles to an¬ 
timonopoles. This is easily understood if we realize that 
the breaking of SU{2)b. x U{1)b-l to U{1)y by (H) and 
(77) is similar to the breaking of the electroweak gauge 
group and, thus, cannot lead to any topologically stable 
monopoles or strings - see Ref. for a more detailed 
argument. It can only lead to the existence of topologi¬ 
cally unstable dumbbell configurations [3l| consisting of 
an open string with a monopole and an antimonopole at 
its two ends. 

The strings, at any given time after their formation, 
can be thought of as random walks with a step of the 
order of the particle horizon [s^ - to describe the evolu¬ 
tion of this string network, we will follow closely this 
reference. They typically connect monopoles to anti¬ 
monopoles, but unstable closed strings of limited size 
may also exist. As shown in Ref. [s^, at all times before 
the entrance of the monopoles into the horizon, there 
is of the order of one string segment per horizon and, 
thus, the ratio of the energy density Ps(t) of the string 
network to the total energy density ptot(i) of the uni¬ 
verse remains practically constant. At cosmic time tu) 
we have approximately one monopole-antimonopole pair 
per horizon volume ~ (47r/3)(2tH)^ connected by an al¬ 
most straight string segment of the size of the horizon 


and energy ^ /is2tH- The energy density Psitn) of the 
strings at cosmic time tn is then ^ 3Gps/2t^. After 
this time, more and more string segments enter the hori¬ 
zon, but the length of each segment remains constant. 
Consequently, the system of string segments behaves like 
pressureless matter and, thus, ps{t) ^ 3Gps/2{tu_t^y^^, 
which implies that the ‘relative string energy density’ 


Psjt) 

p^it) 



(57) 


(p-yit) is the photon energy density) increases with time 
- remember that we are in the radiation dominated era 
of the universe and, thus, ptotit) = p-y{t) = c(T) = 
3/4t^. The strings, finally, decay at cosmic time 

td {TGpg)-^2tB, (58) 


where T ~ 50, by emitting gravitational waves with en¬ 
ergy density Pgw(td) at production given by 


as one can infer from Eq. (IS3. Note that this formula 
also gives the maximal relative string energy density. 

From Eq. (1551) where we substitute the lowest value 
of the number of e-foldings N 2 and take p = 2, we 
find that, for the two numerical examples presented, 
tn ~ 4.76 X 10“^ sec and 1.04 x 10“"'^ sec, respectively. In 
deriving these values, we take g{Tyj = 228.75, which cor¬ 
responds to the spectrum of the minimal SUSY SM, and 
^(Th) = 40.75 and 10.75 in our two numerical examples, 
respectively. The above values of ^(Th) are consistent 
with the effective number of massless degrees of freedom 
at the cosmic temperatures Th corresponding to the val¬ 
ues of the cosmic time tn obtained. We see that the 
strings enter the horizon well before the time of big bang 
nucleosynthesis. Their decay time is ^ 5.97 x 10“^ sec 
and 5.49 sec, in the two cases, as one can find from the 
corresponding dimensionless string tensions 



^^-3.19x10-'^ and 7.57x10"'^. 
2A 


(60) 

This means that the cosmic strings decay around the time 
of nucleosynthesis and, thus, well before recombination 
which takes place at a cosmic time ^ 10^^ sec. As a con¬ 
sequence, they do not affect the CMBR. Their maximal 
relative energy density in the universe is ^ 2.26 x 10“^ 
and 3.48 x 10“^ for our two numerical examples. So the 
cosmic strings remain always subdominant. In particu¬ 
lar, they do not disturb nucleosynthesis at all. 

Had the strings been around until the present time, an 
upper bound would have to be imposed on the dimen¬ 
sionless string tension in order to keep their contribution 
to the CMBR power spectrum at an acceptable level. For 
the Abelian-Higgs field theory model, this bound is found 
to be [s^ 


G/is < 3.2 X 10'^ 


(61) 
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In our first numerical example, the dimensionless string 
tension G^s given in Eq. (150]) almost saturates the upper 
bound in Eq. (IBID , but violates the recent more stringent 
upper bound [S^ 

G/Ts < 3.3 X 10“® (62) 


from pulsar timing arrays, which also holds for strings 
surviving until the present time. Our second numerical 
example violates both the bounds in Eqs. (EB) and (1621) . 
Thus, both our examples are only possible because the 
strings decay sufficiently early. 

The value of the ratio of the energy density of the 
gravitational waves produced by the strings to that of 
the photons at the present cosmic time tg is found from 
Eq. (15^ to be - cf. Ref. - 


Pgvjjtp) 

P'liip) 



(63) 


The present abundance of these gravitational waves is 
then given by 


(to) 


Pgwitp) \ f P'jjto) \ ,2 
P,ito) J \Pc{to) J 


(64) 


where Pc(to) is the present critical energy density of the 
universe and Hq ~ 0.7 the present value of the Hubble 
parameter in units of km sec“^ Mpc“^. We find that, for 
our two numerical examples, r2gw/i^(to) 2.18 x 10“® 
and 3.35 x 10“®, respectively. The frequency /(td) of 
these gravitational waves at production must be ^ 
since the length of the decaying strings is ^ 2tH “ see 
Ref. [13 . The present value of this frequency is then 


where teq is the equidensity time at which matter starts 
dominating the universe. For the two numerical exam¬ 
ples, this frequency turns out to be ^ 1.06 x 10“'* Hz 
and 4.68 x 10“® Hz, respectively. Note that the estimate 
of the frequencies and abundances of the gravity waves 
presented here cannot be made much more accurate, but 
it can certainly be considered good enough for our pur¬ 
poses. 

We see that the predicted frequencies of the gravita¬ 
tional waves are too high to yield any restrictions from 
CMBR considerations [Ig . They are also well above the 
range of frequencies probed by the pulsar timing array 
observations [s^l- So, the most recent stringent bound 
[Isjl from these observations does not apply to our case. 
However, the frequency of the gravitational waves in our 
first numerical example lies marginally within the range 
to be probed by the future space-based laser interferom¬ 
eter gravitational-wave observatories such as the evolved 
laser interferometer space antenna/new gravitational- 
wave observatory (eLISA/NGO) [13, which is expected 
to be able to detect values of Hgwh^(to) as low as 


4 X 10“*°. Our overall conclusion is that the monopole¬ 
string system disappears without causing any trouble, 
but the gravitational waves that it generates may be 
probed by future space-based laser interferometer obser¬ 
vations. 


IV. CONCLUSIONS 

In view of the recent results [l^, [2l| indicating that 
appreciable values of the tensor-to-scalar ratio in the 
CMBR cannot be excluded, we addressed the question 
whether such values can be obtained in SUSY hybrid in¬ 
flation models resulting from particle physics. To this 
end, we have considered a reduced version of the ex¬ 
tended SUSY PS model of Ref. Q, which was initially 
constructed for solving the &-quark mass problem of the 
simplest SUSY PS model with universal boundary con¬ 
ditions. The reason for focusing on this model is that 
it is known to support successful versions of hybrid in¬ 
flation like the standard-smooth one [l^l- This scenario 
is compatible with all the recent data even with a min¬ 
imal Kahler potential, but predicts negligible values of 
the tensor-to-scalar ratio. 

In the context of this particular particle physics model, 
we demonstrated that a two stage hybrid inflationary sce¬ 
nario which can predict values of the tensor-to-scalar ra¬ 
tio of the order of few x 10“^ can be constructed. For 
the values of the parameters considered in this paper, the 
model in the global SUSY limit possesses practically two 
classically flat directions, the trivial and the semi-shifted 
[23 one. The SUGRA corrections to the potential sta¬ 
bilize the trivial flat direction so that it becomes able 
to support a first stage of inflation. All the cosmologi¬ 
cal scales exit the horizon during this inflationary stage 
and our present horizon undergoes a limited number of 
e-foldings. The tensor-to-scalar ratio can acquire appre¬ 
ciable values as a result of sufficiently mild (i.e. non- 
catastrophic), but still appreciable SUGRA corrections 
combined with strong radiative corrections to the infla¬ 
tionary potential, while the value of the scalar spectral 
index remains acceptable. We can obtain values of the 
tensor-to-scalar ratio up to about 0.05. Larger values 
would require unacceptably large running of the scalar 
spectral index. 

The additional number of e-foldings required for solv¬ 
ing the standard problems of hot big bang cosmology 
are generated by a second inflationary stage taking place 
along the semi-shifted path, where U{1)b-l is unbroken. 
This is possible since the semi-shifted direction, being al¬ 
most orthogonal to the trivial path, is not affected by the 
strong radiative corrections on the trivial path and also 
because the SUGRA corrections on the semi-shifted path 
remain mild. 

At the end of the first inflationary stage, the group 
SU{2)b. breaks spontaneously to a U{1) subgroup and, 
thus, magnetic monopoles are formed. The subsequent 
spontaneous breaking of a linear combination of this 
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[/(I) and U{1)b-l at the end of the second inflationary 
stage leads to the production of cosmic strings connect¬ 
ing these monopoles to antimonopoles. At later times, 
the monopoles enter the horizon and the string-monopole 


system decays into gravity waves well before recombina¬ 
tion without leaving any trace in the CMBR. The result¬ 
ing gravity waves, however, may be measurable in the 
future. 
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